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Abstract. We investigate the spectral distortion of the 
cosmic microwave background (CMB) caused by relativis- 
tic plasma. Within the Thomson regime, an exact ana- 
lytic expression for the photon scattering kernel of a mo- 
mentum power-law electron distribution is given, which is 
valid from the non- to the ultra-relativistic regime. The 
decrement in the photon spectrum saturates for electron 
momenta above 3 me c to that of an optically thick ab- 
sorber with the optical depth of the relativistic electrons 
given by the Thomson limit. Thus, the ultra-relativistic 
Sunyaev-Zeldovich (SZ) decrement measures the electron 
number and not the energy content. On the other hand, 
the relativistic SZ increment at higher frequencies depends 
strongly on the spectral shape of the electrons, allowing 
for investigation of relativistic electron populations with 
future instruments. 

We calculate the expected Comptonization due to the 
energy release of radio galaxies, which we estimate to 
be « 3 • 10^^ erg Gpc~'^. We investigate Comptonization 
from (a) the part of the released energy which is thermal- 
ized and (b) the relativistic, remnant radio plasma, which 
may form a second, relativistic phase in the intergalactic 
medium, nearly unobservable for present day instruments 
(presence of so called 'radio ghosts'). We find a thermal 
Comptonization parameter due to (a) of y ~ 10^^ and 
(b) an optical depth of relativistic electrons in old radio 
plasma of Trei < 10"'^. If a substantial fraction of the vol- 
ume of clusters of galaxies is filled with such old radio 
plasma the SZ effect based determination of the Hubble 
constant is biased to lower values, if this is not accounted 
for. Finally, it is shown that a supra-thermal population of 
electrons in the Coma cluster would produce a signature in 
the Wien-tail of the CMB, which is marginally detectable 
with a multifrequency measurement by the Planck satel- 
lites. Such an electron population is expected to exist, 
since its bremsstrahlung would explain Coma's recently 
reported high energy X-ray excess. 
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1. Introduction 

The Sunyaev-Zeldovich (SZ) effect has proven to be an 
important tool for cosmology and the study of clusters 
of galaxies. It measures the total thermal energy content 
of the electron population along a line of sight and does 
not depend on spectral features of the underlying electron 
distributions, as long as they are non-relativistic. For a 
review see Sunyaev & Zeldovich (1980). 

Deviations in the spectral shape of the SZ effect occur 
for higher electron temperatures or energies, increasing 
from the mildly, through the trans- to the ultra-relativistic 
regime. Approximation for relativistic corrections to the 
thermal SZ effect are given in the literature. For reviews 
see Rephaeh (1995a), Fargion & Sahs (1998), Molnar & 
Birkinshaw (1999), Birkinshaw (1999), and Sazonov & 
Sunyaev (2000). 

The thermal SZ effect (thSZ) is therefore a calorime- 
ter for heat releasing processes in the universe, such as 
the gravitational compression of matter during structure 
formation, and the energy output of galaxies. A very vi- 
olent form of energy release of galaxies are the outflows 
of relativistic plasma from active galactic nuclei. This 
plasma fills large volumes, the radio lobes, and thereby 
forms very low density cavities filled with relativistic par- 
ticles and magnetic fields (for evidence for these cavities 
see: Bohringer et al 1995, Carilli & Harris 1996, Carilli & 
Barthel 1996, Clarke et al. 1997; McNamara et al. 2000). 
The replaced inter-galactic medium (IGM) is heated by 
this process (e.g. Kaiser & Alexander 1999). It is one goal 
of this work to estimate the amount of this energy and its 
imprint on the CMB (see also Enfilin et al. 1998b; Yamada 
et al. 1999). 

Also the radio plasma is able to Comptonize the CMB. 
A search for the SZ effect in radio lobes was carried out 
by McKinnon et al. (1991). After being released from the 
radio galaxy, the relativistic electrons with the highest en- 
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Fig. 1. The functions g{x), h{x), —i{x), and above the latter jcr(2;; a, pi,p2) — forpi = 1, 3 (from top to bottom) , 
P2 = 1000, and a — 2. The dashed curve is the relativistic correct, numerically estimated g{x) for a thermal plasma 
with /cTe = 15 keV. 



ergies cool rapidly due to synchrotron and inverse Comp- 
ton (IC) losses, letting the radio lobes fade away after the 
AGN stops its activity. Patches of such remnant, unde- 
tectable radio plasma were named 'radio ghosts' (EnBlin 
1999). However, the less energetic electrons may remain 
relativistic in such an environment for cosmological times. 
Compton scattering off these electrons removes photons 
from the CMB Planck spectrum and scatters them to 
much higher energies, leading to an SZ effect with a char- 
acteristic spectral signature. 

Since radio ghosts are a still speculative ingredient of 
the IGM, we discuss different strategies to reveal their 
presence. We further discuss their possible influence on 
the SZ effect in clusters of galaxies. 

The paper is organized as follows. In Sect. ^ we dis- 
cuss the theory of transrelativistic inverse Compton (IC) 
scattering. In Sect. ^ we examine the possible role radio 
ghosts can have, and estimate their contribution to the 
CMB-Comptonization. In Sect. || we discuss the possible 
influence remnant radio plasma in the ICM can have on 
the determination of the Hubble parameter via the SZ 
effect (Sect. 4.1). Further SZ measurements of cluster of 
galaxies are proposed as a tool to test for non-thermal elec- 
trons in the intra-cluster medium (Sect. ^^ ). Conclusions 
can be found in Sect. ||. 

We assume an Einstein de Sitter (EdS) cosmology and 
Ho = 5O/150 kms~^ Mpc~^. All IC calculations are done 
in the optically thin limit. 



2. Comptonization of the CMB 

2.1. Phenomenology of SZ- Effects 

We distinguish here three processes which change the 
spectral shape of the CMB by inverse Compton scattering: 

— the thermal Sunyaev-Zel'dovich effect (thSZ), 

— the kinetic Sunyaev-Zel'dovich effect (kSZ), 

— up-scattering of photons by relativistic electrons (rSZ). 

The latter happens whenever a CMB photon gets scat- 
tered off a relativistic electron. An electron with gamma 
factor 7 increases the frequency ly of a scattered photon 
{ly') on average by a factor ly'/iy — ^■f^ — i. The pho- 
ton is therefore scattered to much higher energies and is 
effectively removed from the spectral range of the CMB. 
We consider a photon to be removed from the CMB, if its 
energy is increased by about one order of magnitude or 
greater, requiring 7 > 3. Relativistic electrons with lower 
7-factors still produce a substantial energy gain of the 
photon, which also depopulates the Rayleigh-Jeans side 
of the CMB, but these electrons are less effective scatter- 
ers compared to electrons of higher en ergi es. An accurate 
way to treat this case is given in Sect. 2.2. 

Using the usual Kompaneets approximation and ne- 
glecting for the moment those electrons with 7 < 3, the 
spectral distortion produced by the three processes under 
considerations are (see e.g. Rephaeli 1995b) 



5i{x) = 5(a;)2/gas - h{x) f3g 



i{x)T„ 



(1) 
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where the CMB-blackbody spectrum is 



1 



(2) 



with X = hv/kT and iq = 2 {kTcmh)^ / (hc)'^ ■ For compar- 
ison: the peak of the spectrum is at a; = 2.82, and the 
up-coming Planck satellite will measure the CMB signal-^ 
at a; = 1.8, 2.5, 3.8, 6.2, 9.6, and 15 (Bersanelh et al. 1996; J 
Puget 1998; see also Tab. |). The thSZ-distortions, which^ 
have a spectral shape given by ^ ' 



(e^ - 1)2 



1 



1 



(3) 



depend in their strength on the Comptonization parame- 
ter 



dl ric 



(4) 



The line-of-sight integral extends from the last scattering 
surface of the cosmic background radiation dX z — 1100 
to the observer at z = 0. The kSZ-distortions have the 
spectral shape 



h{x) 



(e^ - 1)2 



(5) 



and depend on /3gas, the average linc-of-sight streaming 
velocity of the thermal gas (wgas = /3gas c, /3gas > if gas 
is approaching the observer), and the Thomson optical 
depth 



(6) 



/^gas ^gas — J dl ^T-Cjgas /^gas ■ 

Finally, the optical depth of relativistic electrons 

Trel = CTT I dl Ucrd 



(7) 



2.2. Transrelativistic Thomson Scattering 

We derive here the exact formulae for the photon redis- 
tribution function of the transrelativistic SZ effect in the 
Thomson regime, which can be applied to get the rela- 
tivistically correct thSZ, but also the rSZ for any arbitrary 
isotropic electron distribution. Such a formula has already 
been given in the literature (Rephaeli 1995a), but not in 
the compact form provided here (see Eq. (^l|)). We further 
give the analytically exact photon redistribution function 
for a power-law electron distribution. 

The change in intensity due to IC-scattering of a black- 
body photon distribution at an isotropic distribution of 
electrons in the optical thin limit is given by 



(8) 




S/Smax(p) 



Fig. 2. The frequency redistribution function Ps{s) plot- 
ted in normalized coordinates {Ps{s;p) x Smaxip) over 
s/smax(p)) for different electron momenta p. From the left 
to the right curve is p = 0.1, 0.3, 1,3, yielding Smax(p) = 



0.2,0.6,1.8,3.6 or t^^Up) = 



dp) 



1.2,1.8,5.8,38. 



where t is the Thomson optical depth 

(Tt / dl He 



(9) 



i{x)T is the flux scattered to other frequencies, and 
j (x) T is the flux scattered from other frequencies to 
X — hv/{kT). In Eq. (|l|) we assumed that j{x) i{x) in 
the region of interest {x < 10) for ultra-relativistic elec- 
trons. In the following we drop this approximation. 

Eq. (||) can be re-written to include a y-like parameter 



determines the fraction of the number of photons, which 5i{x) - 
are effectively removed from the CMB. 

where 



y 

kf, 
(kt) 



9{x) V , 



Too 

El 



{j{x) -i{x)) 

Jdl Tie kTf. 
Jdl rio 



dl Ur kTr 



Too c 



(10) 

(11) 
(12) 

(13) 
(14) 



with Pc the electron pressure (see Eq. (A. 6)). The defini- 
tion of the pseudo-temperature kT^ is equal to the ther- 
modynamic temperature in the case of a thermal electron 
distribution (Eq. (p^)). 

If one introduces the photon redistribution function 
P{t) dt, which gives the probability that the frequency of 
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Fig. 3. The CR-frequcncy redistribution function 
Pcr,s{s;a,pi,p2) over s. pi = I, p2 = 10, and a = 
3.0, 2.5, 2.0 from left to right. 

a scattered photon is shifted by a factor i, one can express 
the scattered spectrum as 



dtP{t)i{x/t) , 



(15) 



where P{t) gives the probability that a photon is scat- 
tered to a frequency t times its original frequency. For a 
given electron spectrum /o [p) dp, where p — (3clc is the 
normalized electron momentum, and where Jdp fdp) = 1, 
the photon redistribution function can be written as 



Pit) 



dpUp)p{t-p). 



(16) 



P{t;p) is the redistribution function for a momoener- 
getic electron distribution. P{t; p) can be derived following 
Wright's (1979) kinematical considerations of the IC scat- 
tering in the Thomson regime (ft <C 7e me c^). Scattering 
of a photon from an isotropic photon field has a probabil- 
ity distribution of the angle 9 between photon and electron 
direction in the electron's rest frame given by 



f{^Ji)d^i = [2^t{^-|3,^iY]-'d^l , 



(17) 



where /x = cos 9. The frequency shift from v to v' of the 
photon after scattering into an angle 9' is 



_ j/ _ 1_ 



/3c M 



(18) 



The probability distribution of /i' — cos 9' for a given /i is 
5(m';m) d/i' = ^[1 + AiV" + ^(1 - A^')(l - m")] V (19) 



(Chandrasekhar 1950). The probability for a shift t given 
the electron velocity is 

'dt 



P{t;p) = / dfj.f{fi)g{fi';fi) 



(20) 



fi'(t, (3c, m) is given by ( p^ ) and the range of integration is 
given by the conditions — 1 < /i < 1 and — 1 < /u' < 1. The 
integrand in ( pO| ) is a rational function of and therefore 
analytically integrable: 



P{t;p) = 
3(1 + t) 



3|l-t| 



[1 + (10 + 8p2 + Ap^)t + e] + (21) 
3 + 3p2 + p4 3 _^ 



32pH 



2p 



-(2arcsinh(p) - |ln(i)|) 



The maximal frequency shift is given by 

|ln(t)| < 2arcsinh(p) , (22) 

and therefore P{t;p) = for |ln(t)| > 2arcsinh(p). 
Similar expressions for the photon redistribution func- 
tion using different variables can be found in the lit- 
erature (Rephaeli 1995a; Fargion et al. 1996; Enfilin & 
Biermann 1998; Sazonov & Sunyaev 2000). We verified 
analytically, that P(t;p) has the proper normalization 
{Jdt P{t;p) = 1) and checked numerically that it also re- 
produces the known average energy gain of the scattered 
photons {JdtP{t;p) (1-1) = fp^). 

If one is interested in the large frequency shift {t ^ 
1) due to IC-scattering by ultra-relativistic electrons it is 
convenient to use the logarithmic frequency shift s — hi{t). 
The photon redistribution function then reads 

Psis;p)ds ^ P{e';p)e'' ds, (23) 

and Ps{s;p) = for |s| > Smax(p) = 2arcsinh(p). Psis;p) 
is plotted in Fig ||. 

We are interested in three distinct forms of the mo- 
mentum distribution of the scattering electrons: 

— a thermal electron distribution: 
/3th 



/c,th(p; Ah) = , exp(-/?th VT+F) , (24) 

"2 (/5th) 

with Ki, denoting the modified Bessel function of the 
second kind (Abramowitz & Stegun 1965), introduced 
here for proper normalization, and /3th — fricC^/kTo 
the normalized thermal beta-parameter, 
a power-law distribution between pi and p2 , as a model 
for a cosmic ray (CR) electron population: 

(a - l)p-" 



fc,cr{p;a,pi,p2) 



pI-" 



(25) 



— a thermal distribution with a high-energy power-law 
tail, as a model for a thermal population modified by 
in-situ particle acceleration: 

/o,th&cr(p;/3th,Q;,_pi,p2) = C'c(/3th,a,Pi,P2) x (26) 

/c,th(p;/3th) ; P<Pi 

/c,th(pi; Ah) (p/pi)^" ; Pi<p<P2 
; p2 < p 
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Fig. 4. The CR-frcquency redistribution function in log- 
arithmic units logj^Q(Pcr. s(s; a, pi,p2)) over s. pi = 0.1, 
P2 = 10, and a = 3.0, 2.5, 2.0 from left to right. 

where the normalization parameter Cc(/3th, oi,Pi,P2) is 
determined by / dp/c,th&crb; /3th, a,Pi,P2) = 1- 

The integral in Eq. (p^ ) has to be evaluated numerically 
for the thermal electron distribution. But the analytical 
solution for the power-law case can be expressed as 



Pcr(t;a,Pl,P2) 



-B 



1 



a 



-B 



■3 + a 



-B 



2 ' 
5 -I- a 



l-Q 
Pi 

Q 

^2" 
2 + a 
2~ 

4 -f a 



p^-'^ 16 



(1-f t) X 



+ 



1 - 1 



l+t 



1 + r 



+ 



7 + 3a 

3 + a 
12 + 3Q 

5 + Q 

(2arcsinh(p) - |ln(t)|) 
2^* 



(27) 



+ 



2{5 + a)t 5 + a 3 + a l+a 



V2(t) 



with pi{t) = max(pi, -\/i/2) and P2{t) — max(p2, V^/2). 
Here B^(a, 6) denotes the incomplete beta-function 
(Abramowitz & Stegun 1965), and 



[/(pC=/b2)-/(pi). 



(28) 



Pcr(^; a,?3i,?32) = if |ln(t)| > 2 arcsinh(p2)- This func- 
tion is plotted in Fig. ^ and ^ The resulting change in 
a black-body photon spectrum due to Comptonization by 
CR-electrons is shown in Fig. |l| as the curves labeled with 
j — i. Note that for a given frequency the influence of the 
gain in up-scattered photons of lower initial frequencies, 
jcrix), on the spectral distortions are negligible compared 



Fig. 5. j(x) — i{x) for mono-energetic electron spectra 
with p — 0.3, 1, 3, 10, 30. These curves give the spec- 
tral change per electron of this momentum. Note that the 
absorption feature saturates to —i{x) for p > 10. 



to — i(x), which describes the loss of photons from the 
CMB of intially this frequency, as long as x < 10 and the 
lower cutoff the electron spectrum pi > 3. Therefore the 
rSZ effect can be well approximated as a pure absorption 
process of photons in the CMB-spectral range {x < 10) 
for sufficiently relativistic electrons, justifying Eq. (p. 

2.3. Relativistic SZ- Effect in Emission 

Above the crossover frequency {xcmas > 3.83) the arriving 
up-scattered photons outnumber the disappearing ones 
giving an SZ increment. This has a different spectral shape 
in the relativistic case compared with the classical one. 
The higher the frequency, the stronger is the rSZ effect 
compared to the thSZ effect. This provides for a promis- 
ing opportunity to detect rSZ effect via this characteristic 
signature. Fig. || - illustrate this. A direct application 
of the emission signature of the rSZ effect can be found in 
Sect. O 



3. Radio Galaxies 

3.1. Remnant Radio Plasma 

The presence of radio plasma produces all three SZ ef- 
fects discussed in this paper: The expanding lobes of radio 
galaxies (RGs) and radio-loud quasars compress and may 
shock the ambient IGM, thereby increasing its thSZ effect. 
During the inflation phase of the radio lobes the ambient 
medium is pushed aside, leading to a possible kSZ effect. 
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Fig. 6. The relativistically correct g{x) for mono- 
energetic electron spectra with p = 0.1, 0.3, 1, 3, 10. Also 
the low energy limit g{x) is shown, which is the curve 
with the most pronounced extrema. The curves are pro- 
portional to the spectral change per electron pressure. 



And finally the relativistic electron content of the radio 
plasma produces the rSZ effect. 

The inflation phase of a RG is short compared to cos- 
mological times, thus we do not consider the transient 
kinematic effect. Furthermore, since gas on both sides of 
the radio lobes is accelerated in different directions, giving 
contributions of opposite sign but similar strength which 
should therefore roughly cancel out. 



In Sect. 3.2 we ar- 
gue that the combined kSZ effect of an isotropic ensemble 
of flows is practically indistinguishable from a thSZ ef- 
fect with the same energy content, further justifying our 
simplified treatment. 

Radio lobes become unobservable due to adiabatic, IC- 
and synchrotron energy losses of the electrons on a time- 
scale shorter than the expansion time of the radio lobes 
which leads to pressure equilibrium with the surround- 
ing medium (e.g. Kaiser et al. 2000). However, the energy 
transfered to the IGM by the expansion of the lobes, and 
also the relativistic electron population remain and leave 
their fingerprints in the CMB spectrum via Comptoniza- 
tion. 

The radio plasma and later radio ghosts will expand 
or contract until they reach pressure equilibrium with the 
surrounding medium. The pressure of the radio ghost is 
given by that of the confined relativistic particles and the 
magnetic fields, assumed to be in rough energy equiparti- 
tion. Thus the magnetic field energy densities should be 
of the order of the thermal energy density of the envi- 
ronment. Subsonic turbulence in this environment, which 



Fig. 7. \j{x) — i{x)\ for mono-energetic electron spectra in 
double logarithmic units. The dashed lines show absorp- 
tion, the solid lines emission. The curves are labeled with 
the electron momentum. 



has an energy density below the thermal energy density, 
is therefore not strong enough to overcome the magnetic 
elastic forces of the radio ghost if the magnetic field in the 
ghost is contiguous on large scales. In the case of a mag- 
netic field disjoined on scales smaller than the physical size 
of the ghost and/or the turbulence being sonic or super- 
sonic, which is e.g. expected in giant merger events of clus- 
ter of galaxies, it can 'shred' the ghost into smaller pieces. 
The size of such pieces will be comparable to the eddy 
size of the turbulence. This means, since a typical turbu- 
lent spectrum has less energy density on smaller scales, 
that there is a length-scale below which the turbulence is 
not able to overcome magnetic forces. Turbulent erosion of 
radio ghosts should stop at this length-scale, leaving small- 
scale patches of still unmixed old radio plasma. Similarly, 
fluid instabilities of Kelvin-Helmholtz or Rayleigh- Taylor 
type at the surface of the ghost will not lead to a complete 
mixing of the radio plasma with the surrounding matter. 

The relativistic electrons suffer from adiabatic losses 
during the inflation phase of the radio plasma, after- 
wards only from synchrotron and IC- losses. Coulomb- 
and bremsstrahlung-losses are negligible in a tenuous rel- 
ativistic plasma. The cooling of an ultra-relativistic elec- 
tron with momentum p is governed by 



dp 



ac + abP + asP 



(29) 



(Kardashev 1962), where we ignored adiabatic losses or 
gains due to volume changes. The coefficient gq, a-h, and Og 
for Coulomb, brcmsstrahlung, and synchrotron/IC losses 
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are (Rephaeli 1979, Blumenthal 
3 



Gould 1970) 



ah 



-(It crio In 



0.22 



3a 

TT 

4 (Tt c 



In 2p-- 
^ 3 



3 mo 



(eB + ecmb) , 



(30) 
(31) 
(32) 



where the electron density of the background gas is ne, 
£b = -6^/(8 tt) is the magnetic field and Ecmb the CMB 
photon energy density. The plasma frequency is Up = 
y'47re^no/me and a is the fine- structure constant. In or- 
der to judge the importance of the different terms in Eq. 
( p9| ) we define the cooling times 

- — = 22 Gyr (1 + 0.01 ln(pi/n_5))"' (33) 

= — 8.6 • 10^ Gyr (1 -I- 0.38 In(pi))"^ /n_5, (34) 
ah 



9.8Gyr/(e_nPi), 



(35) 



where pi = p/lO, n_5 = nc/(10~^ cm~^), and e_ii = 
(eb +ecmb)/(10~^"^ ergcm"^). We note, that in the follow- 
ing estimates the radio plasma is assumed to contain rela- 
tivistic electrons only above p = 10. From these timescales 
it becomes obvious that as long as the electron density 
within the radio plasma stays below 10~^ cm~^ only syn- 
chrotron and IC losses have to be taken into account. If 
the radio plasma does not mix on a microscopic scale with 
the ambient medium this can be expected. In the follow- 
ing we assume that this is indeed the case. We note that 
in the case of a cold denser gas component in the radio 
plasma the Coulomb cooling heats the cold electron and 
this heat contributes to the thSZ-effect. The evolution of 
the relativistic electron distribution function in this case 
can be calculated analytically for constant conditions in 
the plasma {B,no, Scmh = const) (Enfilin et al. 1999). If 
Coulomb- and bremsstrahlung-cooling can be neglected 
Eq. (p^) is solved by 



Pit) = 



1 

Po 



3 mc c'^ V 8 vr 



(ecmb) t 



(36) 



where the brackets ((...)) indicate time-averaged quanti- 
ties, and Po is the momentum at t = 0. In order to demon- 
strate that even under disadvantageous conditions elec- 
trons stay relativistic we insert po = 10, {B^) — (20 ^G)^, 
and (Ecmb) = 11 Ecmb.today Corresponding to an time- 
average over the CMB energy density between today and 
z = 2 in an EdS Universe. The final electron momentum 
after t = 10 Gyr coohng is still p > 3. Thus, practically 
all relativistic electrons in old radio plasma stay relativis- 
tic for cosmological times. Their energies change, but this 
is not important for the rSZ effect in the spectral region 



of the CMB: the rSZ decrement depends practically only 
on the number of relativistic electrons, and not on their 
energy. This is because all the up-scattered photons have 
energies far above the CMB range, so that their actual 
energy does not influence measurements within this range 
{x < 10). However, a search for up-scattered photons is a 
promising way to detect the rSZ-effect of very low energy 
relativistic electrons. 



3.2. Detectability of Radio Ghosts 

Beside the possibility to see a rSZ signature, discussed in 
this article, there are a few other possibilities to detect 
radio ghosts. 

Radio ghosts are practically unobservable as long as 
their electron population remains at low energies. But if 
the population is re-accelerated the ghost becomes ra- 
dio luminous again. This can happen when the ghost is 
dragged into a large-scale shock wave, e.g. in a merger 
event of clusters of galaxies or at the accretion shock where 
the matter is falling onto a cluster. The emission region is 
expected to be irregularly shaped, and should exhibit lin- 
ear polarization due to the compression of the magnetic 
fields in the shock. Such regions are indeed observed in 
the periphery of a few clusters of galaxies and are called 
'cluster radio relics' (for reviews, see Jaffe 1992 and Fer- 
etti & Giovannini 1996). Their properties, such as degree 
and direction of polarization, surface luminosity, periph- 
eral position etc., can be understood within such a sce- 
nario (Enfilin et al. 1998a). The observed rarity of the 
cluster radio relic phenomena can be understood in this 
context: the presence of a shock wave, which should be 
quite frequently found in and near clusters of galaxies, is 
not sufficient to produce a cluster radio relic, a radio ghost 
has to be present at the same location, too. 

Another way to detect the presence of radio ghosts 
is via their ability to magnetically scatter ultra-high en- 
ergy (UHE) cosmic rays (CR) (Enfilin 1999). Without such 
scattering the distribution of sky arrival directions of UHE 
CR should trace their source distribution within « 50 
Mpc, due to the limited distance protons above 3 • 10^^ 
eV can travel without strong photo-pion energy losses 
(Greisen 1966, Zatsepin & Kuzmin 1966). If the distri- 
bution of UHE CR sources follows that of the matter in 
the local Universe a non-uniform arrival direction distri- 
bution is expected, contrary to the observations. Medina 
Tanco & Enfilin (in preparation) show that under opti- 
mistic assumptions about their distribution radio ghosts 
can sufficiently isotropize the UHE CR arrival directions. 
Since the UHE CR scattering angle decreases with parti- 
cle energy, this scenario can in principle be tested, as soon 
as sufficient CR data is available. 
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3. 3. The Energy Budget of Radio Ghosts 

The energy a radio galaxy release into a radio ghost and 
its environment during its lifetime is 



h Vgh 



gh 



(37) 



where Vgh is the volume occupied by the ghost, P the 
pressure inside and outside the ghost and 



Ec + + es = (1 + fcp + kB)ec 



(38) 



is the energy density in the remnant radio plasma, split 
into its leptonic (Ec), hadronic (gp) and magnetic {eb) 
part. Unfortunately, the ratios between these energy den- 
sities are still undetermined. The parameters kp = Sp/ee 
and kB = es/se take account of this. 

The radio plasma is expected to follow a relativistic 
equation of state {P = \£gh\ note that isotropically ori- 
ented magnetic fields also follow a relativistic equation of 
state). The fraction of i?gh stored as volume work in the 
ambient medium is therefore 



Ai?th — P Vgti — - Egi, 



(39) 



This thermal energy produces a thSZ effect in addition to 
that expected from the heat in the IGM due to structure 
formation. It seems to be difficult to follow this energy 
in an expanding Universe with ongoing structure forma- 
tion. But it can be argued that this energy is released into 
gravitationally bound structures (clusters and filaments 
of galaxies). Filaments expand due to the Hubble flow, 
but the matter within a given filament flows into the next 
galaxy cluster and gets compressed thereby. This means, 
that on a longer time-scale the extra heat from RG is ex- 
pected to be confined in clusters of galaxies and does not 
suffer from adiabatic losses. 

In fact, part of the energy injected into the environ- 
ment might also be stored gravitationally whenever gas 
is pushed to a higher gravitational potential or trans- 
formed to kinetic energy of gas flows. Gravitational energy 
is converted to kinetic energy, whenever the accelerated 
gas flows into the next potential well. From Eqs. (|l|) and 
(^ one would conclude that the kSZ effect is zero for an 
average of isotropically oriented flows. But Eq. (|l|) is only 
an approximation. A better description of the average kSZ 
effect for an isotropic ensemble of flows can be gained by 
the following argument: the individual electron velocities 
of all flows measured in the CMB-rest frame can be com- 
bined into a single momentum space distribution function, 
which is of course isotropic. This function can be used in 
order to calculate the Comptonization with the formalism 



described in Sect. 2.2, since the latter does not depend on 



the position of an electron along the line of sight in the op- 
tically thin limit. The kinetic energies are non-relativistic, 
so that the spectral changes are well described by the func- 
tion g{x) times the total kinetic energy of the electrons 
(see Fig. ||). Thus the Comptonization is independent of 



the exact spectral shape of the electron spectrum or the 
distribution of flow velocities. Therefore we assume that 
all of the energy, which the inflating radio lobes give to 
their environment, contributes to the thSZ effect. A simi- 
lar argumentation for gaussian, isotropic, random velocity 
flelds was given in Zeldovich et al. (1972). 

We approximate the relativistic electron spectrum by 
a single power-law (Eq. (p^)), so that the energy density 
Sc and the number density rio.cr of the relativistic electron 
of a ghost are related via the average kinetic energy of the 
electrons ((70 — 1) me c^): 

3! — 1 r, 



(70 - 1) TOc C 



a - 



■ rricC pi 



(40) 



The approximation assumes the particles to be ultra- 
relativistic (1 ^ pi < P2) and the distribution to be dom- 
inated by the lower end (pi <^ P2 and a > 2). In the 
case that the distribution has a trans- or non-relativistic 
part the exact formulae for energy density, pressure and 
adiabatic index can be found in Appendix A. 

If the ultra-relativistic approximation is valid and the 
low energy electrons dominate the spectrum, we can write: 

a—1 {1 + kp + kB) rric pi 

All the quantities on the rhs have to be inserted for the 
same stage of the radio plasma evolution in order to give 
a consistent result. The best observational constraints for 
these quantities are of course given for the moment of 
injection, when the radio plasma is visible. 

3.4. CMB- Distortions 

We estimate now the Comptonization due to energy from 
radio galaxies stored in different forms into the IGM. We 
have to calculate the integrals in Eq. (^) and (0) in a cos- 
mological context. We need to know the distribution of 
injected heat and relativistic electrons. Since both quan- 
tities (injected thermal heat in the IGM and number of 
relativistic electrons in the radio plasma) are accumula- 
tive in the sense, that radio galaxies produce them, but 
cooling mechanisms are insufficient to remove them. For 
the thermal energy this is clear, since except for cooling 
flows in the very center of some cluster of galaxies the 
cooling time exceeds the Hubble time. For the number 
density of relativistic electrons this follows from the fact 
that the dominant cooling mechanism in a tenuous rel- 
ativistic plasma are synchrotron- and IC-cooling, which 



become inefficient for low energy electrons (see Sec 3.1). 

The integrands in Eq. (Q) and (0) are proportional to 
the amount of remnant radio plasma at a given epoch (e.g. 
measured by the number density of relativistic electrons 
ne,rei)- Since, as we argued, radio plasma is conserved, its 
electron density is given by 



"•G,rel(0 = / (^^' "-e,rel(i') 



(42) 
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where ne,roi(i) is the source electron density of radio 
plasma at time t in comoving coordinates. The optical 
depth of relativistic plasma can then be written as 



dt 



dt 



Trcl=crTC J — nc.rci(z) J dz' — (1 + z') , (43) 

where t{z) is the time between today and redshift z, and 
-^max is the maximal redshift of injection. We use Zmax = 
4, which is sufficiently low so that cooling effects can be 
neglected. Since no,roi(z), the source density of relativistic 
electrons, is only poorly constrained observationally, we 
use Qjet(i), the source density of energy from RG instead. 
For convenience, we define 

c..h = -^ r"d.^Qjet(.) [dz'^il + z'f, (44) 



rric Jo 



dz' 



and the total amount of energy release through radio 
plasma: 



Cgh 



dz ^Qjct{z) . 
dz 



(45) 



For an ultra-relativistic and steep electron population in- 
side ghosts the approximation Eq. (^) gives 



3 a-2 



A similar expression gives the additional thSZ effect due 
to IGM heating by expanding radio lobes: 

UJgh 



f2 



(47) 



This follows from Eq. (^ if we assume that the IGM 
protons get half of the heat energy. 

The optical depth of radio ghosts can be written as 



a-2 9ygh/j5i 



f (l + fcp + fcs) 



Vgh _ 
10 1 



ks Pi 



(48) 



Here, we inserted a = 2.5. If we assume protons, electrons 
and magnetic fields to have all the same energy density, 
and insert a speculative lower electron cutoff of pi = 10, 
this equation suggests that ygh 3> Tgh. This means that 
the thSZ decrement due to the thermal IGM heating by 
inflating radio lobes always dominates over the rSZ decre- 
ment of the relativistic electron population in these radio 
lobes since \g{x)\ \i{x)\ in the Rayleight- Jeans regime 
(see Fig. |] and Eq. (||) for details). The thSZ increment 
in the Wien-regime decreases exponential and therefore 
much faster than the rSZ increment of a power-law elec- 
tron population. Therefore the rSZ effect of radio plasma 
dominates over its induced thSZ effect at sufficiently high 
frequencies as long as the electron population stays rela- 
tivistic. 

The total jet power of radio galaxies per comoving vol- 
ume can be derived from the radio luminosity function of 
RG with the additional assumption that there is a unique 
relation between radio luminosity and jet power. This is 



not strictly correct, since it is known that for the largest 
part of the observed lifetime of RGs they should exhibit 
some luminosity evolution, even in models with constant 
jet power (for the case of powerful RG of type FRII see e.g. 
Begelman & Cioffi (1989), Falle (1991), Kaiser & Alexan- 
der (1997), Kaiser et al. (1997), Daly (1999)). But this 
evolution is 'modest' (roughly an order of magnitude) and 
its influence should produce some scatter in an empirical 
derived relation, but not a large systematic effect. Sim- 
ilarly we neglect the influence of the density of the gas 
in the radio source environments on its radio luminosity. 
Again we expect this only to increase the scatter about 
the empirical function we will use in the following to re- 
late radio luminosity and jet power. 

EnBlin et al. (1997) have fit a power-law to the jet 
power - radio luminosity relation derived by Rawlings & 
Saunders (1991) for a sample of radio galaxies, which in- 
cludes both, the most powerful radio galaxies of type FRII 
and also the somewhat less luminous FRI objects. Their 
results are based on minimum energy arguments and age 
estimates. The real energy of radio plasma can easily be 
much higher than the minimal energy estimate by some 
factor /power > 1 due to the presence of relativistic pro- 
tons, low energy electrons or deviations from equiparti- 
tion between particle and field energy densities. A rough 
estimate of /power can be derived from observations of ra- 
dio lobes embedded in the intra-cluster medium (ICM) of 
clusters of galaxies. These observations show a discrepancy 
of the thermal ICM-pressure to the pressure in the radio 
plasma following minimal energy arguments of a factor of 
5—10, even if projection effects are taken into account (e.g. 
Feretti et al. 1992). Since also a filling factor smaller than 
unity of the radio plasma in the radio lobes can mimic a 
higher energy density we chose /power = 3 in order to be 
conservative. 

The jet power - radio luminosity correlation a,i v — 2.7 
GHz is 



(?jet(i.) = a, (L,/(WattHz-i h^^)f'^ /, 



power 



(49) 



for which K = 0.82 ± 0.07, logio(a^/ergs-i /igo^) ^ 
45.28 - 26.22 • ± 0.18 (EnBlin et al. 1997). This relation 
allows to integrate the observed radio luminosity function 
n^{L^,t) in order to get the total jet power 



Qjctiz) = J dLy nt,{Lt,, z) (jjet(iz.) • 



(50) 



We use imin = 10^^ Watt Hz^^, which is just above the 
region of the RLF which is dominated by starburst galax- 
ies, and Lniax = 10^^ Watt Hz~^. We adopt different radio- 
luminosity functions parameterized by Dunlop & Peacock 
(1990) and integrate Eq. ( |50| ) in an EdS-cosmology up 
to redshift Zmax = 4. Strictly, the expressions for n^(L^, t) 
used here only apply to those members of the radio source 
population with steep radio spectra (spectral index —0.7 
or less). However, the space density of the flat spectrum 
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RLF 


b. 




Vgh 




10"" erg Gpc"^ 


10"" 


PLE 


0.7 


2.80 


1.51 


PLE 


0.82 


2.19 


1.23 


PLE 


1.0 


1.91 


1.15 


RLF2 


0.7 


4.86 


1.94 


RLF2 


0.82 


3.04 


1.40 


RLF2 


1.0 


2.04 


1.18 



Table 1. Cosmological energy output of radio galax- 
ies for two radio luminosity functions of Dunlop & Pea- 
cock (1990) and expected Comptonization parameter (Eq. 
@). PLE denotes the pure luminosity evolution and 
RLF2 is the second free-form model of Dunlop & Pea- 
cock. The optical depth Tgh of the relativistic electrons in 
ghosts is not well determined since it depends on several 
unknown parameters. Eq. ( p8| ) should be used in order to 
translate ygh to Tgh. If the internal energy of the radio 
plasma is fully thermalized, then T/gh increases by a factor 
of up to 4. 



population, which is thought to consist mainly of FRLtype 
objects, is at least a factor 10 smaller at any given red- 
shift (Dunlop & Peacock 1990). Therefore we can safely 
neglect their contribution to the overall radio luminosity 
function. We further calculate the ghost distribution func- 
tion for biy = 0.7, 0.82, 1 in order to show the dependence 
on the uncertainties. 

Results of the integration of the different RLFs and 
jet power-radio luminosity correlations are given in Tab. 
|l|. The energy input in form of radio plasma is roughly 

Cgh = 3 • 10^6 erg GpC-3 (/power/3). 

Note that the X-ray background, which is believed to 
be dominated by AGN emission, corresponds to an in- 
jection energy of 10^^~^* erg Gpc""^ in comoving coordi- 
nates (Soltan 1982; Chokshi & Turner 1992; Fabian & Iwa- 
sawa 1999; Fabian 1999). Either the X-ray energy losses of 
AGNs exceed the radio plasma release, or if X-ray power 
is comparable to the jetpower then the latter is strongly 
underestimated here and /power ~ 30 would be more real- 
istic. The resulting Comptonization would be comparable 
to the present-day upper limit of y < 1.5 • 10^^ (Fixsen 
et al. 1996). The extra-galactic radio background has an 
energy density of « 5 • 10^^ erg Gpc"'^ (Longair & Sunyaev 
1971), which is several orders of magnitude below the en- 
ergy density in the X-ray background and the expected 
ghost energy density. This indicates that radio emission is 
a very inefficient mechanism to extract energy from radio 
plasma. 

It is interesting to note that the X-ray background 
predicts a mass density in AGN black holes of (1.4. ..2. 2) • 
10^"* Mq Gpc~'^ for a mass-to-light-conversion-efficiency of 
AGNs of 0.1. This is consistent with the central black 
hole mass density derived from the observed black-hole- 
to-galactic-bulge-mass ratio and the observed masses of 
galaxies (Faber et al. 1997). 



If the total energy of the radio lobes would be thermal- 
ized, then no rSZ-effect would result, but the contribution 
to the thSZ effect would increase by a factor of 4, shifting 
it closer to the present upper limit. 

3.5. Discussion 

The expected thermal Comptonization due to heating 
of the environment of RG is of the order of y — 1.4 • 
10~^ (/power/3), if it is assumed that radio ghosts do not 
thermalize their internal energy, otherwise up to a factor 
of 4 higher, /power is the ratio of true energy over equipar- 
tition energy in radio lobes. 

There are two attempts in the literature to estimate 
the thSZ effect due to energy release from RGs. 

Enfilin et al. (1998b) give y < 1.0 • 10"^ ??jet/x-ray For 
this estimate they assume that the jet power is completely 
thermalized, and that the energy budget is the same as the 
X-ray energy release (?7jet /x-ray ~ li thus 10 times higher 
than what we estimate here). It was also assumed that the 
thermal energy suffers from adiabatic cooling due to the 
Hubble flow, but since radio galaxies are located in fila- 
ments and clusters of galaxies only 1-dimensional instead 
of 3-dimentional expansion was taken into account. 

Yamada et al. (1999) use a model for heating and 
cooling of the environment of radio lobes. They esti- 
mate the jet power from the black-hole-to-galaxy mass 
ratio (/bh — 0.002) using a Press-Schechter description of 
galaxy (and black hole) growth and assuming a constant 
fraction (/ig — 0.01) of galaxies above 10^^ M© to be ac- 
tive RGs. Since the total energy release in their description 
should be higher by an order of magnitude, compared to 
what we estimate here, it is not surprising that they get 
y = 5.7-10-5 (/bh/0.002) (/r/0.01), which slightly violates 
the present day upper limit. 

The present day limit on the Comptonization param- 
eter gives already an important constraint on the energy 
budget of RGs, which will be improved in the near future. 
If it would be possible to detect the spectral signature of 
a rSZ effect, a further major step in the investigation of 
radio plasma could be achieved. We do not give detailed 
prediction for the rSZ increment, since different to the rSZ 
decrement at lower frequencies it strongly depends on the 
unknown shape of the electron distribution. But observa- 
tions of the rSZ effect in emission could allow a spectral 
examination of the low-energy end of the relativistic elec- 
tron content of the Universe. 

4. Clusters of Galstxies 

4.1. Embedded Radio Plasma 

Clusters of galaxies are known to be strong sources of the 
thSZ effect. Since the thermal gas causing the thSZ dis- 
tortions can also be observed with X-ray satellites via its 
bremsstrahlung emission it is possible to compare the an- 
gular diameter with the true line-of-sight diameter of the 
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cluster. This gives directly the Hubble parameter, Hg, as- 
suming spherical symmetry at least in a statistical sense. 
It is therefore of principal interest to estimate the possible 
influence of relativistic plasma on the measured Comp- 
tonization. 

Radio plasma was mostly produced by outflows from 
AGN in virialized cosmological structures such as fil- 
aments and clusters of galaxies. Subsequent flows intcg 
deeper gravitational potentials should have transported^ 
a significant fraction of all radio plasma into clusters of 
galaxies. If the radio plasma is still unmixed with the ther- 
mal medium, it resembles a cavity in the X-ray emitting 
ICM gas. We assume in the following that a (for simplicity 
constant) fraction <I>gh of the cluster volume is occupied 
by radio ghosts. The X-ray luminosity of this cluster is 
thus given by 



Lx ^ax dV nl kT,' (1 - $giO - l[ 



1 "c.O 



;i - $gh) . (51) 



Tie is the cluster electron density, rie^o its central value, 
and Zci the characteristic length scale of the cluster (e.g 
the core radius) . Since the cluster temperature kT^ can be 
determined spectroscopically we assume it to be known. 
The thSZ effect 



(jj' 
m, 



^ LnefcTe(l-$gh)-^Zcirio,o(l-$gh) (52) 

■c C J 



in combination with Eq. ( pl| ) allows to measure the typical 
scale of the cluster: 



Ir 



^(1 

Vcl 



(53) 



Since the angular diameter of a cluster should (at least in a 
statistical average) be identical to lc\ the Hubble constant 
can be derived: 



Ho 



Vol 



Lx (1 - $gh) 



(54) 



There are two ways the presence of radio ghosts can af- 
fect the determination of Hq, which shift it in opposite 
directions: First, if a significant fraction of the volume is 
filled with relativistic plasma, then the true Hq is greater 
than the one derived under the assumption that $gh = 0. 
Second, if the measured ?/obs contains a significant con- 
tamination due to absorption by the rSZ effect the true 
yd, and therefore also Hq, is lower. 

The second effect is negligible for the following reason: 
The optical depth of the relativistic electrons of the ghosts 
can be written as 



6 (jjgh Vcl 



(1 + fcs + kp) (7e - 1) (1 - $gh) ' 



(55) 



if we assume pressure equilibrium between ghosts and 
thermal plasma, and use Eqs. (^), (|5^), and |A.^ . Op- 
timistic values for the unknown parameters of the old 
radio plasma (1 -I- fc^ -I- fcp « 3, (7c — 1) ~ 10) give 



Fig. 8. The models 0-3 for the electron spectrum in Coma 

/c,cl(p)- 



Tgh ~ 0.2 $gh j/ci- The rSZ decrement is roughly — Tgh i(x), 
or smaller. The observed, rSZ-contaminated y-parameter 
is therefore 



— t (*-^) 

2/obs = Vcl H ^ Tgh 



j/ci(l + O.O6 0gh). 



(56) 



In the last approximation we assumed that the measure- 
ment is taken around x — h v /kTcmh ~ 2, where the thSZ 
absorption is at its maximum and —i{x)/g{x) <, 0.3 (see 
Fig. |l|) . Inserting ?/obs instead of yd into Eq. (|^ gives an 
systematic error of ~ 0.1 $gh, which is an order of magni- 
tude smaller than the error due to the factor 1 — ^gh in the 
denominator of that equation. The error due to the rSZ 
effect could be larger, if the measurement is done closer to 
the crossover frequency at x = 3.83. However, this value 
of X is inconvenient for detecting the thSZ effect. Also if 
(7c — 1) w 1 the distortion would be significant, but this 
is extremely speculative. 

We conclude that the presence of a significant fraction 
of the volume of the ICM filled by radio ghosts biases the 
determination of the Hubble constant to lower values, if 
this effect is not accounted for. 



4-2. Non-Thermal Electrons in the ICM 

In some clusters of galaxies a relativistic ICM population 
of electrons is visible due to their synchrotron emission, 
which produces the radio halos of clusters of galaxies. Also 
the recently detected extreme ultraviolet (EUV) excess 
(Lieu et al. 1996) and the high energy X-ray (HEX) excess 
(Fusco-Femiano et al. 1998, 1999) in the Coma cluster 
indicate the presence of non-thermal electrons. 

Our knowledge about the slope of the electrons spec- 
trum in the ICM is limited (see Enfilin & Biermann 1998 
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Fig. 9. log^o \5Iu/ (Jy arcmin^^)! over i^/GHz for the cen- 
ter of the Coma cluster. The solid lines are models 0-3 
and the dashed lines are the difference of models 1-3 to 
model 0. The dots give the frequencies and sensitivities 
of single beams of the Planck satellite (Bersanelli et al. 
1996; Puget 1998; see also Tab. ||). Since the angular di- 
ameter of Coma (core radius =10') is much larger than 
the beam of Planck (FWHM < 10') a marginal detection 
(see Appendix B) of a possible non-thermal electron pop- 
ulation producing the HEX excess by bremsstrahlung can 
be expected. Balloon experiments or other high-altitude 
observatories (at the planned ALMA site or in Antarc- 
tica) will be much more sensitive. Note that the curves on 
the Ihs represent absorption, whereas the ones on the rhs 
represent emission. 



for an attempt to compile the electron spectrum in the 
Coma cluster). The number density and therefore the op- 
tical depth of the higher energy electrons is too small to 
detect the rSZ decrement produced by them. See Birkin- 
shaw (1999) for a discussion of the rSZ decrement from 
radio emitting electrons and McKinnon et al. (1990) for 
an attempt to detect this decrement in radio galaxies. The 
rSZ effect emission from these electrons may already have 
been observed in the EUV (Hwang 1997, Enfihn & Bier- 
mann 1998, Sarazin & Lieu 1998, Bowyer & Berghofer 

1998) or in the HEX excess (Fusco-Femiano et al. 1998, 

1999) . The latter would require a high number of electrons 
in the range of 3 — 5 GeV, which would over-produce radio 
emission if the magnetic fields are as high as Faraday rota- 
tion measurements indicate (2— 10 /LtG; Crusius-Waetzel et 
al. 1990; Kim et al. 1991; Feretti et al. 1995, 1999, Clarke 
et al. in preparation). Only if fields are as weak as 0.16 fxG 
consistency between the expected and observed radio flux 
is established (Fusco-Femiano et al. 1998, 1999). As a pos- 



sible solution of this apparent contradiction of magnetic 
field estimates Enfilin et al. (1999) and Sarazin & Kemp- 
ner (2000) proposed that the HEX excess could alterna- 
tively be produced by bremsstrahlung of a non-thermal 
high energy tail of the thermal electron distribution. Such 
a tail could exist due to in-situ particle acceleration pow- 
ered by ICM turbulence (Enfilin et al. 1999; Dogiel 1999; 
Blasi 2000). 

Detection of the IC emission of such a non-thermal tail 
could confirm the bremsstrahlung-scenario. In the follow- 
ing we estimate the expected spectral distortions for the 
Coma cluster of galaxies. 

The electron density of the cluster is assumed to follow 
a beta-profile: 



nc{r) = nc,o [l + {r/rc\Y 



(57) 



where rd — 400 kpc h^Q^ is the core radius, — 0.8 is the 
beta-parameter, and n^fi = 2.89 • 10"'^ cm~^ is the central 
thermal electron density (Briel et al. 1992). We ignore the 
possible complication due to embedded old radio plasma 
discussed in Sect. LI for simplicity. 

The optical depth of a line of sight passing the cluster 
center at a distance R is then 



Tcl(i?) 



1 



[1 + {R/rcx) 



(58) 



giving a central optical depth of rci(O) = 5.95 -10 ^ . The 
CMB-distortions are 



5Iy = io 5i{x) = io {j{x) - i{x))) . 



(59) 



The electron spectrum is given for a thermal spectrum by 
Eq. ( p^ ) and for a modified thermal spectrum by Eq. ( p6| ) . 
The non-thermal electrons are treated as an additional 
population, increasing slightly the total electron number, 
but leaving the number densities of electrons with p < pi 
unchanged. We estimate the distortions for a pure ther- 
mal spectrum (model 0) and for three modified thermal 
spectra, which are able to explain the HEX excess of 
Coma by bremsstrahlung (see Fig. 6 in Enfilin et al. 1999): 

model 1: a = 5.3, pi = 0.45, i?kin4 = 50keV; 

model 2: a = 2.9, pi = 0.51, i;kin4 = 62.5 keV; 

model 3: a = 1.6, pi = 0.58, Ei^i^^i = 80keV; 
P2 = 10 in all three models. The electron spectra are 
shown in Fig. || and the resulting CMB distortions at 
the center of the Coma cluster in Fig. ^. The differences 
between the thermal and the modified thermal models 
are considerable and change as a function of frequency, 
so that future multichannel CMB telescopes such as 
high-altitude ground based observatories (e.g. at the 
ALMA site or in Antarctica), balloon-experiments, or 
marginally even the Planck satellite (see Appendix B) 
can discriminate between these models. A detection of 
this non-thermal SZ effect would prove in-situ acceler- 
ation processes acting in the ICM of Coma. It would 
confirm the bremsstrahlung origin of the HEX excess and 
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therefore solve the discrepancy between Faraday and IC 
based magnetic field estimates in Coma. 

Since the SZ effect does not depend on distance (as 
long as the instrument beam resolves the cluster core) 
many other clusters can be investigated for the presence 
of a non-thermal SZ effects with future CMB experiments. 



5. Conclusion 

We investigated the transrelativistic Thomson scattering 
of photons on a isotropic electron distribution in the opti- 
cally thin limit using the photon redistribution kernel (Eq. 
(^)). We derived for the first time an analytic formula for 
the scattering by a population of electrons with a power 
law momentum distribution, Eq. (p6|). We demonstrated 
that relativistic electron populations produce a decrement 
in the cosmic microwave background, similar to that an 
absorber with an optical depth of the Thomson depth 
of the relativistic electrons would produce. Our formal- 
ism can be generalized to the optical finite case using the 
techniques described in Rephaeli (f995a,b) or Molnar & 
Birkinshaw (f999). 

We applied this theory to radio galaxies and clusters 
of galaxies: 

Although a single radio galaxy produces a negligible 
SZ decrement, the combined effect of several radio galax- 
ies and their remnants might produce a detectable sig- 
nal. In order to show this we estimated the total cos- 
mological release in jet power of radio galaxies using the 
observed radio luminosity function converted by an em- 
pirical jet power-radio luminosity relation. It is roughly 
3- 10*5'^ erg Gpc~^ (/powor/3), where /power gives the poorly 
constrained ratio between true energy content of radio 
lobes and the minimum energy estimate. If completely 
thcrmalized this energy would lead to a Comptonization 
parameter ofy~6-10~^, close to the present day obser- 
vational limit of y < 1.5 • 10~^. We argued that roughly 
3/4 of the released energy remains as a relativistic plasma, 
which rapidly becomes unobservable after the activity of 
the galaxy stopped. Thus the Comptonization due to heat- 
ing by radio plasma is expected to be y « 1.5 • 10^^. 

Patches of old plasma, called 'radio ghosts', are ex- 
pected to survive turbulent erosion unmixed with the am- 
bient medium, so that they are able to retain a low en- 
ergy relativistic electron population. The optical depth 
for Thomson scattering by radio ghosts is Tgh ~ 10"^ for 
very optimistic assumptions about the low energy cutoff 
of the fresh electrons during injection. Otherwise Tgh is 
much lower . 

If radio ghosts occupy a significant volume in clusters 
of galaxies, they affect the determination of the Hubble 
constant via SZ and X-ray measurements. The geomet- 
ric effect of the cavities formed by ghosts in the intra- 
cluster gas overwhelms the SZ decrement expected due to 
up-scattering of CMB photons by the ghosts' relativistic 



electron populations. SZ based Hq determinations have to 
be corrected to higher values due to this effect. 

Finally, we demonstrated that future CMB telescopes 
such as the Planck satellite seemed to be useful tools to 
measure supra-thermal electron populations in clusters of 
galaxies. These are expected in the case of turbulent in- 
situ particle acceleration, and supported by the recently 
detected high energy X-ray excess in the Coma and Abell 
2556 cluster. Such a detection would be crucial for our un- 
derstanding of particle acceleration processes in clusters. 

Acknowledgements. We thank Rashid Sunyaev, Matthias 
Bartelmann, and Mark Birkinshaw, the referee, for comments 
on the manuscript. 



Appendix A: Energy Density and Pressure of a 
Trans-Relativistic Power-Law 
Electron Distribution 

If an electron population, with number density richer, has 
a momentum spectrum which is a single power-law (Eq. 
(p6|)) the kinetic energy density Sc is 



oo 



"o.cr / dp fc{p) ( Vl+P^ - 1) mo C 
/o 

„2 



^^c.cr m^ C 



2 i+p^ 



a — 2 3 — a 



1 b'-"]|i 

pi ^c,cr C 



a - 2 
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(A.3) 
(A.4) 
(A.5) 



Here, we have used the short-hand notation defined in Eq. 
(p8|). The first approximation assumes the particles to be 
ultra-relativistic (1 <C pi < P2), and the second that the 
distribution is dominated by the lower end [pi <^ p2 and 
a > 2). For an ultra-relativistic population of electrons 
(and similar for protons) the pressure is Po = \^c, but 
this is not correct for a transrelativistic population. There 



Pc = 
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The ultra-relativistic approximation was again applied, as- 
suming 1 <^ pi < p2- If this is not justified, correct results 
can easily be obtained from Eq. (A_^) and ( A. 7 ). 



14 



T. A. Enfilin & C. R. Kaiser: Comptonization of the CMB by Relativistic Plasma 



Channel k 




1 


2 


3 


4 


5 


6 


Frequency 


GHz 


100 


143 


217 


353 


545 


857 


Frequency Xk 




1.76 


2.52 


3.82 


6.21 


9.59 


15.1 


Beam FWHMfc 


arcmin 


10.7 


8.0 


5.5 


5.0 


5.0 


5.0 


Sensitivity AT/T 


^.K/K 


1.7 


2.0 


4.3 


14.4 


147 


6670 


Sensitivity Ai(a::fc) 


10"^ 


4.09 


7.66 


21.0 


43.2 


85.2 


97.6 



Table 2. The observing frequencies and sensitivities of the Planck satellite (Puget 1998). Note that the sensitivity 
AT/T refers to fluctuations in the radiation temperature, rather than to the system temperature of the receivers. 



Appendix B: Detectability of the 

Cluster-rSZ-EfFect with 
Multifrequency Experiments 

In order to demonstrate that the predicted non-thermal 
CMB distortions in the clusters like the Coma cluster 



(Sect. 4.2) are detectable with multifrequency experi- 
ments we estimate the expected sensitivity of the Planck 
satelitte. The spectral distortions at the location of the ob- 
served cluster are produced by temperature fluctuations of 
the CMB on the scale of the cluster {ST/Tk 10"^... 10"^), 
the thSZ and kSZ effect of the ICM gas, and the rSZ ef- 
fect of the non-thermal electron distribution. The resulting 
distortion are therefore 



with 



(B.l) 



Tl = Tth, fl{x) = jth{x) ~i{x), 

T3^ST/T + 0^,,, h{x)^h{x). 

Here, jici{x) are the spectral distorions produced by 
the electrons in the non-thermal tail with optical thick- 
ness Troi and spectrum given by /e.th&crb; /3th, a,pi,p2) - 
fc,th{p', Pth) (with proper normalization). 

The observing frequencies and sensitivities of the chan- 
nels (fc) of the Planck satelitte are taken from Puget (1998) 
and given in Tab. |. We assume the beams to be Gaus- 
sian, and therefore to have the areas Ak — 1.13 FWHM^.. 
The core region of the Coma cluster is oleoma ~ ti" 7'corc ~ 
300 arcmin^. Thus, for each frequency Nk — Acoma/Ak 
independent beams probe the CMB. Using a analysis 
allows to dissentangle the different spectral distortions if 
sufficient multichannel sensitivity is given. The accuracy 
in the determination of tj is given by Atj = 
C = A-^ and 



Cjj, where 



. ,r fj{xk)fi{xk) 



(B.2) 



(Press et al. 1992). Assuming Tth — 5.95 • 10^^ and Tici 
and jrci according to the models 1,2, and 3 (see Sect. 4.2) 
we find that Ti-ci/Atioi — 0.61, 1.15, 0.96 correspondingly. 
This demonstrates that the Planck mission is expected 



to give a marginal (1-sigma) detection of such a supra- 
thermal electron population in the Coma cluster. Since 
some other clusters are also expected to contain supra- 
thermal electrons - e.g. Abell 2256 revealed very recently 
a similar HEX excess (Fusco-Femiano 2000) - a combined 
signal from several (HEX excess selected, or merger and 
post-merger) clusters might be detectable with statistical 
significance by Planck. This demonstrates the ability of 
future sensitive multichannel CMB telescopes ~ e.g. dedi- 
cated baloon experiments - to detect supra-thermal elec- 
tron populations in clusters of galaxies. 
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